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We show that instanton calculations in QCD become theoretically well defined in the 
gluon saturation environment which suppresses large size instantons. The effective cutoff 
scale is determined by the inverse of the saturation scale. We concentrate on two most 
important cases: the small- x tail of a gluon distribution of a high energy hadron or a large 
nucleus and the central rapidity region in a high energy hadronic or heavy ion collision. 
In the saturation regime the gluon density in a single large ultrarelativistic nucleus is high 
and gluonic fields are given by the classical solutions of the equations of motion. We show 
that these strong classical fields do not affect the density of instantons in the nuclear wave 
function compared to the instanton density in the vacuum. A classical solution with non- 
trivial topological charge is found for the gluon field of a single nucleus at the lowest order 
in the instanton perturbation theory. In the case of ultrarelativistic heavy ion collisions a 
strong classical gluonic field is produced in the central rapidity region. We demonstrate that 
this field introduces a suppression factor of exp {— c p 4 Qt / [8 a 2 s N c (Q s t) 2 ] } in the instanton 
size distribution, where Q 3 is the saturation scale of both (identical) nuclei, r is the proper 
time and c ~ 1 is the gluon liberation coefficient. This factor suggests that gluonic saturation 
effects at the early stages of nuclear collisions regulate the instanton size distribution in the 
infrared region and make the instanton density finite by suppressing large size instantons. 



I. INTRODUCTION 



At very high energies corresponding to extremely small values of Bjorken x variable the density of partons 
in the hadronic or nuclear light cone wave function can get very large leading to the effect of saturation 
of gluon and quark distributions p]. The large density of partons in the transverse plane produces strong 
gluonic fields ~ l/g which leads to a number of interesting non-linear phenomena. The transition to the 
saturation regime in the (x,Q 2 ) plane is characterized by the saturation scale Q s (x). It has been shown in 
||| H that for a large nucleus the saturation scale grows with the nuclear atomic number as Q 2 ~ A 1 / 3 , and 
thus for a heavy ion it can get large (Q 2 3> &q CD ) making the strong coupling constant small a s (Q s ) <C 1. 
This allowed McLerran and Venugopalan [Q to argue that gluonic fields in the saturation region are given 
by the solution of the classical Yang-Mills equations of motion in the presence of the source given by the 
ultrarelativistic nucleus on the light cone. The solution of the classical Yang-Mills equations for a large 
nucleus (non-Abclian Weizsacker- Williams field) has been found in pM. It has been shown |^|,^J|] that 
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saturation regularizes the power-law divergence of the classical unintegrated gluon distribution function in 
the infrared region. The unintegrated gluon distribution at the two-gluon level is proportional to 
and diverges for small k±. The classical field resums all multiple rescattering effects, which corresponds 
to summing up powers of Q s /k±. The resulting unintegrated gluon distribution given by the non-Abelian 
Weizsacker- Williams field of a large nucleus is proportional to lnQ s /k± in the small transverse momentum 
region Thus the infrared singularity becomes only logarithmic and integrable. As we go towards 

smaller values of x quantum corrections to the classical field become important. The summation of leading 
logarithmic (In l/x) corrections to the classical multiple rescattering picture led to a non- linear evolution 
equation for the total scattering cross section of a color dipole on a hadron or nucleus which was derived 
independently by one of the authors in Q using Mueller's dipole model of [|l0| and by Balitsky in 
employing the high energy effective lagrangian technique. Similar equation has emerged recently out of 
renormalization group approach |l^ , |l3|] . The solution of this non- linear evolution equation should give the 
behavior of the hadronic cross sections at very high energies and specify the dependence of Q s on x. The 
solution has been found by approximate analytical methods |) 14 and by numerical simulations |l5|,[lF 



yielding us with the saturation scale which is a growing function of energy, Q s ~ l/x s , where S is close to 



the value of the BFKL pomeron intercept |17|. Thus at very high energies the saturation scale can get large 
even for a hadron making the small- x tail of the gluon distribution in hadrons similar to the small- x tail of 
the gluon distribution in large nuclei. 

With a large fraction of high energy scattering data being attributed to some non-perturbative QCD 
phenomena, instantons [[lq-pifl might play an important role in high energy scattering processes |22] . There 



have been developed several techniques of calculating the instanton-mediated scattering amplitudes including 



the effective instanton lagrangian approach 23 24J and the valley method |2jJ. Instanton-induced effects 
could be important to deep inelastic scattering processes and other high energy processes as was shown in 
pp| p8| . Several years ago there has been developed a vigorous activity studying the possibility of baryon 
number violating effects due to instantons in electroweak theory |2!| . Recently the authors suggested that 
instanton-induced interactions could play a key role in the dynamics of soft pomerons in hadronic scattering 
cross sections and in multi-particle production (3l) (see also f^). 

The role of instantons in heavy ion collisions has not been studied in such detail. It has been recently 
suggested |3^| that the instanton-induced particle production can account for a significant fraction of the 
total charged particle multiplicity in the heavy ion collisions at RHIC . Since it has been argued that the 
small- a; gluons are most important for the mid-rapidity particle production in heavy ion collisions |35|-|3^| it 
is natural to investigate whether instanton-induced effects are important for these small- x gluons, which is 
our main goal here. 

The qualitative picture of what one might expect to happen to instantons in the saturation environment 
is as follows. Let us consider a gluon propagating through the strong background gluonic field created either 
by a single large nucleus or in a collision of two nuclei. The gluon would undergo multiple rescatterings in 
the field which would have the physical effect of generating some non-zero effective mass for the gluon, which 
would be roughly proportional to Q s (see, for instance, |39| ). Equivalently we can say that due to multiple 
rescatterings the gluons are screened on the transverse distances inversely proportional to the effective mass, 
i.e., on the distances of the order of 1/Q S . Here we can draw an analogy to the case of QCD at high finite 
temperature T. There the Debye mass of the gluons in quark-gluon plasma mo ~ gT introduces screening 
on the distances of the order oil/ ma ~ 1 / gT which in turn suppresses large size instantons by the factor 

ofEr} 



-const' p 2 T 2 



(1) 



In these approximate estimates we are not keeping track of the factors of g, Nc and Nf, putting them to be 
some constants. Thus arguing that the saturation scale plays the role of Debye screening mass for small- a; 
gluons mr> ~ Q s @ we might expect that the instanton size distribution in the saturation background 
would also obtain a suppression factor of 

g — const p 2 Q 2 

for large instantons. This would imply that instantons with sizes p > po ~ 1/Q S are exponentially suppressed 
and only small instantons with sizes p < pa can contribute to the scattering processes. This would justify 
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the perturbative calculation we are about to perform below, since the strong coupling constant for small size 
instantons is small a s (p < 1/Q S ) <C 1. 

Suppression of large size instantons could be anticipated based on the similar arguments related to the 
philosophy introduced in [Q : the typical size of the nuclear color charge fluctuations in the transverse plane is 
of the order of 1/Q S - Thus instantons of larger size would just be washed out by the color charge fluctuations 
in the ultrarelativistic nucleus. 

An analogy can also be drawn with QCD at finite density and zero temperature. There the instantons are 
again suppressed by an exponential factor of Q 

e ~const p 2 p 2 ^ 

The suppression is due to the following qualitative picture: large size instantons tend to produce quarks 
with momenta k ~ 1/p < fi and this process is suppressed at large /i due to Pauli blocking since all the 
quark energy levels with momenta k < fj, are occupied []. Instantons in the saturation environment produce 
gluons, which are of course bosons and do not have any Pauli blocking effects. The role of Pauli blocking 
for gluons is played by the distribution of gluons in the saturation wave function of a single nucleus which 
is very similar (at the qualitative level) to the Fermi-Dirac distribution. In the case of nuclear collision 
the transverse momentum spectrum of produced gluons also exhibits saturation behavior and levels off in 
the infrared resembling Fermi-Dirac distribution pq). Due to multiple rescatterings the gluons with lower 
transverse momenta get pushed towards higher transverse momenta of the order of k± ~ Q s . Thus there are 
very few gluons left with small k± and the large size instanton-like solutions which would tend to produce 
those should be suppressed. Associating the width of the nuclear gluon distribution with the chemical 
potential of the Fermi-Dirac distribution Q s ~ /i we would again arrive at Eq. (^) . 

The intuitive arguments presented above suggest that the interaction with gluons in the saturation envi- 
ronment will exponentially suppress the instantons with large size (larger than 1/Q s (x)). Such suppression 
would result in an infrared safe and, therefore, well defined theoretically instanton calculus in the saturation 
medium. To show that this is indeed the case is the main objective of this paper. 

The paper is organized as follows. In Sect. II we present the general techniques that we are going to 
use. We will briefly review the solution of the problem of small instantons in a slowly varying background 
field using the effective instanton lagrangian approach |23|,|24|]. We show that the background field affects 
the instanton size distribution only if it has a non-zero value of (or G^G^) at the space-time point of 
interest [||]2j]. 

We will proceed in Sect. Ill by addressing the issue of instanton size distribution in the small- x tail of 
the gluon distribution in a single nucleus or hadron. Throughout most of this paper we will consider the 
classical picture of saturation. We are interested in the small- x tail of the distribution with x small enough 
for multiple rescattering effects to be become important |^,^ || . For this we require the coherence length of a 
small- x gluon to be comparable to nuclear size l c ~ l/2mx > 2R, which yields us with x < 1/AmR p|, |43p3| ] . 
At the same time the value of x should not be too small since we do not want the quantum corrections to 
start playing an important role. Thus we want a s Inl/x < 1 which leads to x > exp(— l/a s ). In Sect. IIIA 
we demonstrate that the classical gluon field of a single hadron or nucleus has zero value of G°^ v and does 
not affect the instanton size distribution, introducing no enhancement or suppression of instanton effects. 
We proceed in Sect. IIIB by calculating the one loop leading logarithmic (In X/x) correction to the classical 
field of the nucleus. For this correction to become important we have to relax the a s InX/x < 1 condition. 
We observe that even after inclusion of this correction the field strength squared G a £ v of the gluon field still 
remains zero and does not affect instanton distribution. We relate this observation to the representation of 
quantum evolution as a series of classical emissions as advocated in JTo| , p~2|] a nd based on that argue that our 
result is true to all orders in evolution resumming powers of a s InX/x |[[ll). 



The classical field of a large nucleus found in |||7]] was a non-topological solution, in the sense that it did 
not create a transition from a region of one topological charge to a region of a different topological charge. In 
this paper we would like to address the issue whether there exists a topological solution of the classical Yang- 
Mills equations of motion in the presence of external source given by the large nucleus. In Sect. Ill we show 



1 Yu. K. would like to thank Larry Yaffe for providing this argument. 
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that this solution, if exists, is just as probable as the usual "non-topological" solution of |6|,|7[. Unfortunately 
we were unable to construct an exact instanton-like solution for this problem. Instead in Sect. IIIB we 
are going to consider a QCD instanton perturbatively interacting with the classical field of the nucleus HQ 
and producing a combined classical field of the instanton-nucleus configuration which would now connect 
regions of different topological charge. This field squared contributes to the one loop correction to G a2 , and, 
if viewed as a field of the nucleus probed by a point-like instanton is not a classical field anymore. However 
the field itself is a classical field of the instanton-nucleus configuration. This classicality is due to the fact 
that we will be performing the calculation resumming all powers of the parameter a 2 A 1 / 3 (or, equivalently, 
Q^/fcjJ and neglecting higher order corrections in a s to it MgL The equivalence between the classical field 
techniques and this resummation has been discussed in p[]6| |8|,|35| p8| and references therein. 

In Sect. IV we will consider the case of hadronic or nuclear collisions. As has been argued in |35|] the 
dominant particle production mechanism for central rapidity gluons could be due to the strong classical 
gluon field produced by the colliding hadrons or nuclei. The field is again characterized by the saturation 
scale of the colliding nuclei Q s . This gluonic field has been calculated at the lowest order in perturbation 
theory in p5| , ^6| and analyzed numerically in [|37|. Recent progress in calculating the produced particle 
spectrum due to the classical field was made in |38| . The classical field produced by colliding nuclei has a 
non-vanishing value of G a2 v and therefore influences the instanton size distribution. By a direct calculation 
using the lowest order classical field J|5|,|36| we demonstrate in Sect. IV that the effect of this classical fields 
would be to suppress large size instantons by a suppression factor of 



8o2JV c (Q,t„) 

with the proper time of the instanton position To = ^/2xo+xo- and c the gluon liberation coefficient. £0 
is the space-time point where we measure the instanton density. Eq. (|J) shows that at large proper time 
(r — > 00), long after the collision, the lagrangian density G a £ v at each particular point gets small and large 
instantons will not be suppressed anymore. The amount of suppression depends on the gluon liberation 
coefficient, which was estimated numerically to be c = 1.29 ± 0.09 |37) and analytically to be c « 2 In 2 |38|l , 
while RHIC data |54| suggests c = 1.23 ± 0.20 p5[ . We will also propose that the instanton suppression of 
Eq. (W) at extremely high energies might lead to suppression and, therefore, unitarization of the soft pomeron 
of [[nj. We will conclude the paper by estimating that the saturation effects reduce the instanton density 
in the central rapidity region at RHIC by three orders of magnitude compared to the instanton density in 
vacuum as extracted from lattice data WA- 



II. INSTANTONS IN BACKGROUND FIELD 



The problem of small instantons in a slowly varying background field was first addressed in [£3|,Q and was 

resolved by introducing the effective instanton lagrangian L^hx). The complete field of a single instanton 
solution could be reconstructed by perturbatively resumming the powers of the effective instanton lagrangian 
which corresponds to perturbation theory in powers of the instanton size parameter p 2 . In our case here 
the background field arises due to the strong source current J^. The current will be due to a single nucleus 
(Sect. Ill) and two colliding nuclei (Sect. IV). Perturbative resummation of powers of the source current 
term translates itself into resummation of the powers of the classical field parameter a 2 . A 1 / 3 Thus the 

problem of instantons in the background classical gluon field is described by the effective action in Minkowski 
space 

Seff = J d 4 x [~G^{x)G%{x) + Li ff (x) + L T eff (x) + J^A^x)^ . (5) 
In Eq. (||) the point-like instanton vertices are given by the instanton-induced effective lagrangian |23|j24| , |4^] 



Li ff (x ) = J dpn (p) dR exp f-^j- p 2 ^ v R aa ' G^ v (x ] 
in which no(p) is the instanton size distribution function in vacuum given by p 



(G) 
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0.466 e - 1679JV <= 1 / 2tt ^ 2N ° 



no{p)= (N c -mN c -2)\7\^i) (7) 

where b = (11/3) N c — (2/3) Nf. In Eq. (^|) Tjff„ v 1S the 't Hooft symbol in Minkowski space denned in terms 
of the usual 't Hooft symbol in euclidean space by Eq] 



n M 



Va^u, /-*, v = 1,2,3 

(8) 

^o4v) A* = 0, v = 1,2,3. 



xo in Eq. (j^) is the position of the instanton and R aa is the matrix of rotations in the color space with 
dR denoting the averaging over instanton color orientations. To obtain the effective lagrangian for anti- 
instantons from Eq. one has to change rj^ 1 into r^^f in it (r)%L, = (5j2u„)*). 

The classical current of a single large ultrarelativistic nucleus in the effective action of Eq. (0) is given in 



McLerran-Venugopalan model by 



Jl = S,^8{x-)p a {xx), (9) 



where p a (x^) is the two-dimensional color charge density of the nucleus. In the calculations of diagrams 
below we will be using the explicit model of the nucleus as consisting of independent nucleons which is 
justified at high energies and is equivalent to the description of the nucleus in terms of the current of Eq. (^) 
HU . In the case of two colliding nuclei the source current will be 

Jl = S p+ 6(x-)pUx ± ) + V 6(x+) fi(x x ), (10) 

with pf(x±) and p%(x±) the color charge densities of the colliding nuclei. 

The action of Eq. (^) solves the problem of topologically non-trivial classical fields of nuclei at the con- 
ceptual level. It allows one to construct a perturbative series in the powers of p 2 and a 2 A 1 / 3 which, with 
proper regularization of singularities, should sum up to yield us the classical field of a nucleus or nuclei with 
non-zero topological charge. In Sect. IIIB we will construct an example of such field in the single instanton 
sector at the lowest order in these parameters. 

Following |2^,^| we may write the first correction in p to the instanton size distribution generated by the 
background field of the nucleus or nuclei as 

n sat (p) = (A\L[ ff \A) = 



/ 3 4 \ 

= «o(p) f 1 + Qs( ^f_ 1} (A\ Gl v {x ) - G^(xo)G^(x ) \A) + higher orders in/ J (11) 

where G®„(xo) is the field strength of the slowly varying background field taken at the position of the 
instanton and G a pv = (l/2)e ilvp(T G a p(y is the dual field strength. (A\ . . . \ A) denotes averaging in the nuclear 
wave function(s) |^,^|. Eq. ( p"l| ) could be obtained by expanding the effective lagrangian Lf^ in the powers of 
p 2 and averaging the resulting terms in the nuclear wave function and over instanton orientations [p3| , p4|]4§| l . 
As was shown in Q for the classical field of a large nucleus at the leading powers in A the higher order 
correlators factorize (Gaussian averaging) yielding the property which was assumed for the vacuum fields in 

@ 



(2k 1)!! [(^| rj™ v R aa ' G°'„ fj^ R w G% \A) 



(12) 



Using the color neutrality of a nucleus we write 
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sa'b' 



(A\ G^G b af} \A) = j^j( A \ G%G a af} \A). 



(13) 



Employing Eq. (|l3|) in Eq. ( |l2] ) together with the orthogonality of the color rotation matrices R aa R ba = 5 ab 
and performing some simple algebra of 't Hooft symbols one can show that the effect of higher order terms 
in p 4 in Eq. (|ll]) is just to exponentiate the lowest order term leading to (cf. |23| , p4|| ) 



n S at(p) = n (p) exp 



TT 3 p 4 



a s (N 2 - 1) 



(A\ G a {xq) G a (xq) - G a (xo)G a (xq) \A) 



(14) 



Thus in order to find the effect of a particular background field on the instanton size distribution all one has 
to do is to calculate the matrix element in the power of the exponent in Eq. (|l4|). The effective instanton 
lagrangian approach is strictly valid only for instantons of the size much smaller than the typical variation of 
the external field [g4| though it usually works for larger instantons too. In the case of saturation background 
there is only one dimensional parameter characterizing the external field and it is the saturation scale Q s . 
Thus one might argue that the matrix element of G a2 is proportional to Q 4 and therefore Eq. (|lj) resums 
all powers of the parameter p 4 Q 4 and is strictly speaking valid only when 



P 4 Qi « i- 



(15) 



This is the approximation that we are going to employ throughout the paper: the instantons we consider are 
much smaller than l/Q s - This assumption also keeps a s (p) small allowing us to use perturbation theory. 

Eq. (48) includes only the dipole interaction term between the instanton and the background field. Higher 
order multipole interactions in general should also be included. However, since these interactions would 



correspond to subleading in A correction to factorization of Eq. (12) the multipole terms would also be 
subleading in A. They would have less powers of A 1 / 3 per power of p 4 than the leading (dipole) term and 
could be neglected. 



III. SINGLE NUCLEUS CASE 
A. Lowest Order Diagrams 

Let us start constructing the modified instanton size distribution in the external field of a single nucleus 
along with the field of the instanton- nucleus configuration using the action of Eq. (||). We will consider a 
single instanton interacting with the classical field of a large ultrarelativistic nucleus. At the lowest orders in 
p 2 and g the field of the I-nucleus configuration would be given simply by the sum of the fields of the instanton 
and the nucleus, similarly to the sum ansatz for the I— I configurations |5l}| . Since we are interested in the 
effect of the nuclear gluon field on the instanton size distribution it is more convenient to start analyzing 
diagrams for the interaction of the instanton with the field directly, without first deriving the field. This is 
equivalent to calculating the matrix element in the exponent of Eq. (fl4|). We will also concentrate for now on 
the case of a single rescattering in the nucleus, which parametrically corresponds to the case of a 2 A 1 / 3 < 1. 
We will generalize our results to the a 2 A 1 / 3 ~ 1 case after we obtain the lowest order expression. 

All the graphs that we are going to analyze will be calculated in d^A^ = light cone gauge. The lowest 
order diagram which might contribute to the action of the I-nucleus configuration is the one gluon exchange 
diagram depicted in Fig. [l]A There an instanton interacts with one nucleon in the nucleus which is moving 
ultrarelativistically with the large momentum in the light cone "+" direction. To obtain the average value 
of the action we have to average this diagram (as well as all others) over the nuclear wave function Q^-j^] ■ 
The procedure includes averaging over all possible positions of nucleons in nuclei and quarks in nucleons, as 
well as averaging in the color spaces of each particular nucleon, which makes sure that the nucleus is in the 
color singlet state both before and after the interaction 0] . Averaging in the color space on the nucleon in 
figure Fig. |l|A gives zero for this diagram since TrT a = 0. 

Similar color arguments can rule out a number of diagrams with extra gluons involved where only one 
gluon interacts with the nucleon. Diagram in Fig. [I]B has two gluons interacting with the nucleon and can 
not be proved to vanish by simple color algebra. Adding the diagram with the gluon lines crossed to it and 
performing color averaging one can see that the diagram is proportional to 



G 



S((p - I - l'f) R aa ' Va>-» I, R bb ' rtf_ v I S ab (uip - I ~ 07+ ^'} P l)2 fr- l+uip) 

+ u(p-l- 07+ Jl^I^J +ftfr)) (16) 

where p is the momentum of the quark line in the nucleon below having a large p + component. All the 
other light cone components of the gluon's momenta in the diagrams we consider are much smaller than that 
l+,l' + <C p+ [|| . Using the orthogonality of the color rotation matrices R aa R ab — 8 ah and evaluating the 
products of 7-matrices in the eikonal approximation one can show that Eq. ( |lt3| ) is proportional to 

^i,*r£-ui'u5(i- + L) ( T J —-T^r-) ~ = °- ( 17 ) 

p \t_ — te t_ + le J 

Thus we have shown that the contribution of the diagram in Fig. to the action of the I-Nucleus system 
is also zero. 




ABC 

FIG. 1. Lowest order diagrams which might contribute to the instanton action in the field of a nucleus. Only one 
nucleon participates in the interaction here. After a simple calculation (see text) one can see that the contributions 
of these diagrams are zero. 

We have to point out that due to the non-Abelian structure of the effective instanton lagrangian of Eq. (^) 
there is another type of an instanton vertex with two gluons fusing into an instanton at the same point 
in coordinate space ]30[| , which yields us with another diagram which might contribute to the instanton- 
nucleus interactions as shown in Fig. |l|C. However the vertex is antisymmetric in the color indices of the 
gluons connecting to it. Due to color averaging in the nucleon below the diagram in Fig. [l|C is proportional 
to f abc 5 ab — and does not contribute to the action. 

The reason why all the diagrams in Fig. |l| are zero could be understood in terms of the effective instanton 



lagrangian result given by Eq. (14). At the leading (classical) order the value of the matrix elements of the 
local operators (xo)G* (xo) and G? I/U (xo)G? IiU (xq) in the nuclear wave function could be obtained by just 
substituting the field strength of the classical non-Abelian Weizsacker- Williams field of a nucleus in it. As 
could be seen from the exact expression for this field given in Q the only non-vanishing components of its 
field strength tensor are G^ , so that {G% Wa f = G^ u Wa G^ Wa = 0. Therefore from Eq. @ it follows 
that there is no correction to the instanton density if one calculates the matrix elements of G 2 and GG only 
at the classical level. However, there exists a possibility that higher order in a s corrections might yield a 
non-zero value for these matrix elements. To check this we will now go one step beyond and calculate what 
appears like a one loop correction to this classical result. 



B. Next-to-Leading Order Diagrams 



In this subsection we are going to analyze two-loop diagrams contributing to the action of the instanton- 
nucleus configuration. These diagrams are of the order of p in the instanton size and involve only one 
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interacting nucleon, that is only one power of a^A 1 / 3 < 1 . First let us point out some two-loop diagrams 
which vanish or can be shown to be suppressed after a simple calculation. Those include the diagrams 
where only one gluon interacts with the nucleus and the diagrams with the vertices of the type introduced 
in Fig. [j]C. For instance the graph in Fig. ||A vanishes for the same reason as the graph in Fig. [l]C: it is 
proportional to f abc S ab = 0. All the diagrams generated by including all other possible connections of the 
extra gluon in Fig. to the gluon and quark lines vanish for similar reasons of simple color algebra. The 
case when the gluon connects to the instanton via the instanton-two gluon vertex is a little different. To 
estimate the diagram in Fig. |^B we have to make use of the fact that the nucleus is moving very fast in the 
light cone "+" direction. Then the gluon propagators for the exchanged gluons give g v _ in covariant 
gauge. The contribution of the graph in Fig. 0B is therefore proportional to 



n n n M R dd ' f acd ' n M R ee ' f bce ' X ab 
9p- 9v- Vdfip n J Veup n 



-3N c g_ 



0. 



(18) 



The gluon labeled p in Fig. gB could also connect the instanton-two gluon vertex to the other two gluon 
lines. These kind of diagrams involve one instanton-single gluon vertex and one instanton-double gluon 
vertex, as shown in Fig. ||. The diagrams in Fig. || can also be safely neglected for the following reasons. The 
graph in Fig. is zero since in covariant gauge its contribution is proportional to Tj^^, g^- g v - ~ g = 0. 




+ 
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FIG. 2. Some of the two-loop graphs involving the instanton-two gluon vertex that vanish. 
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ABC 

FIG. 3. A set of the two-loop graphs involving one instanton-double gluon vertex and one instanton-single gluon 
vertex that vanish. 

To prove that the contribution of the diagram in Fig. |jp3 is zero one has to add to it the diagram where 
the lines corresponding to I and V gluons connect to the nucleon in the opposite order (crossed diagram), 
which similarly to how in was shown in Eq. ( |l7)) for the graph of Fig. would give us two delta functions 
5{l-) S(l'_). Thus one can show that diagram of Fig. is proportional to 
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V^A(k-2l)^- + (-2k + l)- 9liV + {l + k)^] v A a [ _ p l< p 5{l_)5{l'_) 



= -3rJ^_J vV ^_ p l' p S(l.)5(l'_) = 3l-l'_5(l,)6{l'_) = 0. 



(19) 



Finally one can show that in the eikonal approximation employed here the contribution of the graph in 
Fig. |^C is independent of l + . Therefore the integral over l + in that diagram yields us with || 



dl+ 
~2^ 



,(x - 



-) - 



6{x - 



) 



(20) 



where xq is the position of the instanton and Z- is the (frozen) light cone coordinate of the interacting quark 
in the nucleon. Eq. (^5|) requires that the nuclcon interacting with the instanton at point xq should have the 
same light cone coordinate as the instanton. But this would not allow us to sum over all nucleons situated at 
different light cone coordinates at the given impact parameter to obtain the enhancement of this diagram by 
the powers of atomic number A 1 / 3 [pPj35|]. In other words the diagram of Fig. 0C would not be resumming 
powers of the parameter a 2 A 1 / 3 , it would have the factor of a 2 without any A / 3 enhancement. Thus it is 
zero at the leading order in A 1/>3 . 
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FIG. 4. Definition of the effective triple gluon vertex in d^A^ = covariant gauge. 



To start analyzing the diagrams with two instanton-single gluon vertices let us define an effective triple 
gluon vertex as shown in Fig. || All the diagrams on the left of Fig. [| are taken in covariant gauge with the 
gluons' light cone momenta being much smaller than the quarks light cone momentum. The graphs in Fig. ^ 
correspond to the first non-trivial perturbative contribution to the combined field of the instanton-nucleus 
configuration. An analogy could be drawn with the Lipatov vertex being the first non-trivial contribution 
to the field of two colliding nuclei fl3"5|,ppS|] . After a simple calculation one can see that the lowest order field 
of the I-nucleus system is 



A;(k) = gf abc (T c ) 



R 



bh' 



4ttV 



d 4 i 



o —il-(x —z) 



(27r)<J(Z_) 



(fe-0c 



1 



{k-lf + ie k 2 + 



ic 



VV-a 



k- 



(21) 



where we Fourier transform over the momentum I flowing between the nucleon at z to the instanton at xq. In 
Eq. ( pl| ) (T c ) is the SU (3) matrix in the color space of the interacting nucleon and I denotes the transverse 
component of momentum I. 

To calculate the contribution to the action at the order p 4 of the field in Eq. ( ^l| ) we have to consider the 
diagram depicted in Fig. ||A. After averaging in the color space of the nucleon below || the graph in Fig. |^A 
yields us with 



1 n 2 

2 y(4^ p ) 



d 4 k dH d 4 l 



4// 



-U'(xq—z) —U'-{xq — z) 



(2tt) 4 (2tt) 4 (2tt) 4 
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X 



{k — l) a (k + l) a < (1 . _ M , . _ M , _ M , _ M 1 



(k-iy + ie {k + l'Y-ie \l 2 1 '' b - av " - Jlt , " a " rKm ' ' J 1 ''D-a»A.- A _ 

1 [- % M a , (fc + 2i% - VbL' (k - l%'9»- + V^a' 2A_] + 5/*-TT~ J ( 22 ) 



x 



where the field on the right hand side of Fig. |5|A is complex conjugate to the field on the left hand side. The 
first factor of 1/2 in Eq. ( |22| ) is the symmetry factor while the second factor of 1/2 arises after averaging over 
colors in the nucleon. Simplifying the expression in Eq. ( |2^ ) by employing the algebra of 't Hooft symbols 
p9[ we obtain 

,s A . ^j^^L^-u •<„-„-«■,..-., Mj(i _ )(W j_zi_ 



[-13fc 2 + 14*-/- 14fc • Z' + 24^ • V - U 2 - 5{' 2 ]. (23) 



(fc - 2 + ie (k + I') 2 - ie ff 




B 



FIG. 5. Diagrams with two instanton-single gluon vertices that cancel each other at the leading order in A. 

To average over the nuclear wave function one has to average over the positions of all the quarks in the 
nucleons and nucleons in the nucleus and sum over all the nucleons and quarks in them ||. This could be 
summarized by the following operation 



E E 

nucleons val. quarks 



d 2 z dz- 



(24) 



which has to be applied to Eq. (23). For simplicity we consider a cylindrical nucleus with cross sectional 
area S±_ and the longitudinal extent of 26_ in the infinite momentum frame. Applying Eq. (24) to Eq. ( p3| ) 
yields for a quarkonium nucleus of p,p(§ 



i (S A ) = -lA^p^g 2 — 



, , 2 2A f d 4 k d 2 l dl+dl' + 



S± J (2tt) 4 (2tt) 2 (2tt) 2 k 2 [(k~l) 2 + ie][(k + l') 2 - ie] (l 2 f 



x [-13fc 2 + 14fe • I - Uk ■ I 1 + 14{ 2 ]. 
Performing similar calculations for the diagram in Fig. 0B we obtain 



(25) 
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_ 4 4 2 2^1 /• ri 4 fc rf 2 Z e -^o-^»- 
H B; ~ P ff Sj_ 7 (2tt)4 (2tt)2 (2tt) 2 [(fc-0 2 + ie][(fc + Z') 2 -^] (I 2 ) 2 ' 

Adding up the contributions of Eqs. ( pq ) and ( |2q ) we write 

4 4 2 2A /■ d A k d 2 l dl + dl'^ e -U+x --il' + x - k 2 



i (Sa+b) = -«47r p g 



S± J (2tt) 4 (2tt) 2 (2tt) 2 fc 2 [(fc-/) 2 + ie][(fc + /') 2 - ie] (I 2 ) 2 

x (-7)[(k-l) 2 + (k + l') 2 }. (27) 

If one wants to perform the 1+ and l' + integrals in Eq. ( p7f ) one has to pick up the poles given by the [k — I) 2 
and [k + I') 2 denominators. However each of the terms in the square brackets at the end of Eq. ( |27j ) cancels 
one of these denominators making the expression in Eq. (j2^) zero for any non-zero xq- . Thus the diagrams 
in Fig. H cancel each other and give zero at the leading powers in A. 

That way we have shown that even at one loop order the field of a large nucleus does not affect the 
instanton distribution 

nttiP) = MP)- (28) 

The physical reason behind this result of our calculation is the following. We have seen in Sect. IIIA that 
for a purely classical field the only non-zero component of the field strength tensor is G^J^ a ^,0], which 
leads to G 2 — GG — and due to Eq. (|lj) the field does not affect the instanton size distribution. The 
diagrams we have analyzed above, especially the graphs in Fig. correspond to including one rung of the 
QCD evolution in energy. (We were interested in the real part of the diagrams and did not get the factor 
of hxl/x from them.) Recently, in [ ]Io| , |l2"[ a picture of the evolution in Inl/x has been developed which 
represents the quantum evolution as a series of classical emissions. At each step of the evolution the existing 
partons act as classical sources of color charge and emit a gluon. The gluon then gets incorporated into the 
source and acts as a color charge emitting gluons in the subsequent steps of the evolution^. This tells us 
that if we look at the gluon field generated by any number of evolution steps it would always look like a 
classical field emitted off some complicated (evolved) source moving in the "+" direction and therefore has 



only G 1 ^^ 11 non-zero component of the field strength tensor. Thus from Eq. (14) we derive that this field 
can not change the instanton density. Multiple rescatterings, i.e., higher powers of a^A 1 / 3 would not help 
to obtain a non-zero effect since G 2 = GG = for the full non-Abelian Weizsacker- Williams field of a large 
nucleus which includes the effects of all the multiple rescatterings in it 



i 




FIG. 6. Last diagram which along with the graphs in Fig. ^ contributes to the topologically non-trivial classical 
field of a large nucleus at the lowest order. 



2 This picture is based on the achieved long ago understanding that the appearance of ln(l/:r) for each emitted gluon 
reflects the fact that this gluon lives much longer light cone time than all gluons emitted after it but much shorter 
time than gluons emitted before it |47|. 
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Since the field of a single nucleus does not suppress instantons it would be interesting to construct a 
solution of the Yang-Mills equations of motion with the nuclear source carrying some non-zero topological 
charge. The effective action of Eq. (||) allows us to do that perturbatively. The lowest order field is just a 
direct sum of the fields of the instanton and the nucleus. At the lowest non-trivial order in p 2 and g almost 
all the field is given by the diagrams contributing to the effective vertex of Fig. || and is written down in 
Eq. . Another diagram which contributes to the field at this lowest order is shown in Fig. [| and involves 
and instanton-double gluon vertex. 

Adding the contribution of the diagram from Fig. |^ to the contribution of graphs from Fig. ^| we obtain 
the following expression for the total field of the instanton-nucleus configuration at the lowest order in 
momentum space 

A%{h) = gf abc (T c )R bb '47r 2 p 2 J «.-«•(»"-*) (2n)6(l.) 



(k-l)„ 
{k-l) 2 +ie 



-jzVvvc [(k - 2J)„s„- + (/ + k) u g„. 



2k. 



9^1 



ri M 
Vb'- 



— n M 
,2 



(29) 



To summarize the results of this section we once again point out that the saturation effects fail to either 
enhance or suppress instantons in a single hadron or nucleus (Eq. (|2q)), which makes possible the existence 
of a classical solution for the gluon field of the nucleus carrying a non-zero topological charge, a lowest order 
expression for which is given by Eq. (|2^ ) for the case of unit topological charge. Now that we have demon- 
strated the formalism for calculating the effects of a background classical field on the instanton distribution 
we continue by considering the case of hadronic or nuclear collisions. 



IV. HADRONIC AND NUCLEAR COLLISIONS 

It was shown in |35l |38fl that the gluon production in the central rapidity region of a heavy ion collision at 
very high energies is dominated by the classical gluonic field produced by two colliding nuclei. The colliding 
nuclei could be visualized as ensembles of point color charges moving without any deflection through each 
other along the light cone |35) , similarly to the case of a single nucleus in McLerran- Venugopalan model Q . 
As the nuclei pass through each other the color charges in each nucleus get rotated in color space by the field 
of another nucleus, which leads to a bremsstrahlung emission of gluons off these charges |?6). The gluons 
then subsequently multiply rescatter in the background fields of both nuclei before being produced |B3| . 

The produced classical field is boost invariant and fills the whole region between the nuclei moving apart 



after the collision [B5 38 . The field is produced at all the impact parameters where the collision happens. 



Thus for central collisions the classical field covers the whole nuclear cross sectional area in the transverse 
direction S±_ . In the longitudinal direction the typical gluon production time is of the order of 2k + jk 2 . 



Since for classically produced gluons |fe| ~ Q s and fc+ ~ Q s [|38| , |4lU53| the typical longitudinal time is of 
the order of 1/Q S - This translates into the typical thickness in the z-direction of the region in which the 
classical fields are produced being of the order of I ~ 1/Q S - As was argues by Baier et al in ^] even though 
the classical field, as a leading order in a s term, would still exist at later times the thermalization effects 
would become important there significantly modifying the distribution of produced gluons. In our effective 
lagrangian approach we consider relatively small instantons with sizes p <C I ~ l/Q s (see Eq. (|l5|)) which 
might easily appear in the spatial region described above at times of the order r ~ 1/Q S - For such small 
instantons it is reasonable to consider the problem of instantons in the classical field produced in a nuclear 
collision before thermalization effects became important. 

To estimate the effect of this classical gluon field produced in mid-rapidity on instantons one has to 
calculate the field strength tensor of this field and substitute it into Eq. (|1J) averaging over the wave 
functions of both colliding nuclei. An analytical expression for the classical field exists only at the lowest 
non-trivial order j3^j3(|. It is possible to calculate the multiplicity distribution of the produced particles 
analytically to all orders in the background field without deriving an explicit expression for the field f38fl . 
Numerical results exist for the field and multiplicity distribution |37|. We begin by calculating the field 

strength tensor at the lowest order corresponding to the case when a^A\' 3 < 1 and a 2 A^ 3 < 1, i.e., the weak 
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field case. A\ and A2 are the atomic numbers of the colliding nuclei. In this case only one nucleon from 
each nucleus contributes to the classical field. We will later generalize our result to the strong field case of 



ajA{ 



1/3 



1 and a 2 A 



1/3 



1. 



+ 



+ 




y x i?Xi _ 



FIG. 7. Instanton in the background gluon field of colliding nuclei, which is taken at the lowest non-trivial order 
in g. Thick dots represent Lipatov vertices. 

The lowest order classical gluon field of two colliding nucleons i and j having coordinates 2^ , Xi- and y , 
Uj + correspondingly in covariant gauge d ■ A — is given by 

e 



(2tt) 4 k 2 + iek Q J * (2tt) 



fj 9" rabc /rpb\ /rficN i[*+aU- +k- y J+ - k- y ~q- (x^ - y . )] G^(k,q) 



(30) 



where C^{k,q) is the Lipatov vertex W\ 



(fc-g) 2 
k+ + ze 



k_,2q-k 



(31) 



In Eq. (31) the four- vector is shown in terms of its components in the (+,—,_!_) form. (Tf) and (TJ) 
in Eq. (30) are matrices in the color spaces of the colliding nucleons |],|3(|. Regularization of the gluon 
propagator k 2 in Eq. (|30j) corresponds to using retarded Green function and is used to insure casuality of 
the classical field: it can be non-zero only in the forward light cone |$(| . 

Let us start by calculating the averaged value of the square of the field strength tensor ([G^fso)] 2 ) to 
be used in Eq. (jlj) at the lowest order in a s . One can easily see that at the lowest order in the coupling g 
only the Abelian part of would contribute. We are interested only in the ([G^(xo)] 2 ) in the forward 
light cone in the central rapidity region. Thus we do not need to include the effects of the lowest order 
fields of each of the nuclei which give non-zero G" only on the light cone (at x + = and/or x_ = 0). 
The diagram that we need to calculate is shown in Fig. [?] and corresponds to the lowest order gluon field 
(order g 3 ) interacting with the instanton at the lowest order in p 2 (order p 4 ). Blobs in Fig. denote Lipatov 
vertices. 

The contribution of the diagram in Fig. [7] is 



(g;„( Xq )g;„(x )) lo = ([d^M-dvAKxo)} 2 ) 



(32) 



with xq the position of the instanton. Substituting the field of Eq. 
colors and positions of nucleons using Eq. (G4f) we obtain 



Oh into Eq. (B2h and averaging over 



LO 



AiA 2 

(2tt) 2 S 1± S 2 _ 



2C f 



d A k dk' dk' e -ik + x _-ik_x 0+ -ik' + x _-ik'_x 0+ 



(2tt) 4 (2tt) 2 



(k 2 +iek )(k' 2 + iek' Q ) 



c ^72 

J [ q 2(k~q) 2 } 2 ^ Cu{k ^ ] ~ kv K C ^ k '> -l) - K C,{k', -£)] 



(33) 
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where throughout Eq. ( |33| ) we imply that U = —k. In arriving at Eq. (|33|) we for simplicity assumed that the 
colliding nuclei are cylinders in the z direction with cross sectional areas of Si± and S 2 ± correspondingly. 
As was shown in ||36|| in the integration over the light cone components of k and k' only the poles in the 
propagators k 2 and k' 2 contribute. If one picks up the poles in the Lipatov vertices (see Eq. (j3l])) the 
resulting contribution would be non-zero only on the light cone (36). However here we are interested in 
the central rapidity region only and that contribution would not be important to us. Thus only the gluon 



propagator poles would contribute in Eq. (33). Anticipating this we may estimate the value of the terms in 
the square brackets at the end of Eq. ( |33"| ) employing the k 2 = and k' 2 — conditions. After some lengthy 
algebra one ends up with 

[k„ C v {k,q) - k v C„{k,q)] [*£ C v {k', -q) - k' v C^k' , -q)]\ k 2 =k ,2 =0 = 

= ^ [-2{k ■ qf + k 2 q 2 + 2(k ■ q)q 2 (q 2 ) 2 ] . (34) 



Substituting Eq. (|3J) into Eq. (|3J) yields 

9 



LO (2tt) 2 Si^S 21 _ F J (2^) 4 k 2 [q>{k- q) 2 ] 2 



{GIMGIM) LQ = ^^^2C F 



dk+dk-dk' + dk'_ e -ik+x --ik-x 0+ -ik' + x --ik'_x 0+ 

k ■ k . (35) 



(2tt) 2 (fc 2 + iek ){k' 2 + iek' ) 



One can readily check that 



[-q 2 (k- q) 2 + 2k 2 q 2 -2(k-q) 2 } _ 

[q 2 (k-q) 2 } 2 ~ fc 2 ' 



d TTn ^12 = — = ~72- ( 36 ) 



To obtain Eq. ( p6|) it is easier to first integrate over the angles between fc and q after which the integration 
over \q\ becomes trivial. Employing Eq. ( |36| ) in Eq. ( |35| ) we arrive at 

( G ^( x o)G a ^(x )) = - V 1 2 2Cf8tt I 



'lo ( 27r ) 2 S 1± S 2± J (2^) 4 (fc 2 ) 2 

dk+dk-dk'.dk' e -ik+*o--ik-x 0+ -ik' + x --ik'_x 0+ 

(2n) 2 (fc 2 + iefc )(fc' 2 + iek' Q ) ' 1 ' 

Evaluating the integral over the light cone components of fc and fc' in Eq. ( |37|) we note that the integral over 
k±_ in Eq. ( |37| ) is dominated by small transverse momenta fcj_. Taking the fcj_ — ► limit of the longitudinal 
integral we obtain 



/ 



dk + dk-dk' + dk'_ e -ik+xo—ik-x +-ik' + x —ik'_x + ^ 1 

W? (fc 2 + iefc )(fc' 2 + ^) k-k^-z, (38) 



where tq = ^/2xo+xo- is the proper time. Here we made use of the fact that we are interested in instantons 
in the forward light cone xq+ > and xo- > 0. The approximation used in obtaining Eq. ([58]) is equivalent 
to taking the limit of early proper time, which more formally means tqQ s <C 1. This approximation is not 
crucial for our approach and is made only to simplify the calculations. 
Employing Eq. ( |38| ) in Eq. ( |37| ) yields 

3 A X A 2 „ !6 1 f d 2 k 



(GU*o)GZM) LO = -a» C F -,-,J^. (39) 
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Before we discuss the issue of regularization of the divergent integral over transverse momentum in Eq. ( p9| ) 
let us first evaluate the other correlator in the exponent of Eq. ( |l4| ) . The lowest order field of Eq. (|3^) yields 
after averaging over nucleons similarly to Eq. ( p3|) 

/ - \ g 6 AlA2 [ d i k dk ' +dk '_ e -ik+x --ik-x 0+ -ik' + x --ik'_x + 

[G^(x )G^(x )) lo - — ——2C F 



lo (2tt) 2 S 1± S 2 ± J (2tt) 4 (2tt) 2 (fc 2 + iefc )(fc' 2 +iek' ) 



J [g2~7fc — g) 2 ] 2 ^ e f ll/ P (T ^v(^3 2) C a {k , — £j 



(40) 



where again fc' = — fc. After some simple algebra one obtains 

£ t ,u P ak fl C v (k,q)k' p C a (k',-q)\ k 2 =k ,2 =Q = 8 e^ v k^ cfc q ■ (fc - q), (41) 

where we have again anticipated that the integration over the longitudinal momenta has to pic k up the poles 
of the gluon propagators giving fc 2 = fc' 2 = similarly to the way we used it in obtaining Eq. (|34j). Eq. (|d]) 
together with Eq. ( f40| ) gives 

(g^(xo)G«M) lo ~ J [q2{ f_ q q)2]2 ^^q^q-(k-q) = 0. (42) 



The fact that the integral over q± in Eq. (42) is zero could be seen by changing variables q — > fc — q which 
would demonstrate that the integral is equal to its negative. We conclude therefore that at this lowest 
non-trivial order in the coupling constant 



GlM G U x o)) LO = 0- (43) 

We think that ([43| ) is a general result required by the symmetry of the problem. Indeed, since the topological 
charge d45| ) determines the net helicity of the system, and on the average the net helicity generated in the 
collision is equal to zero, the expectation value of the topological charge should vanish. This does not, 
however, imply a zero dispersion in the topological charge distribution. The event-by-event fluctuations of 
topological charge can exist, and would induce parity-odd correlations in the multiparticle production |p^| . 

The only non-zero correlator for the classical gluon field produced in a hadronic or nuclear collision is 
thus (G° J/ (a;o)G l ° 1/ (a;o))- The correlator in Eq. ( p9[ ) has an infrared-divergent integral in it. The divergence 
is similar to the infrared singularity present in multiplicity distribution of the produced gluons given by the 
lowest order perturbative diagram in Fig. ^ (without the instanton). As was argued in p5|- |38| , |53t multiple 
rescatterings would regularize the integral by effectively inserting an infrared cutoff in the fcj_ integral in 
Eq. ([39]) which would be proportional to the saturation scale Q s . Unfortunately we do not know exactly the 
coefficient of proportionality between this effective cutoff and Q s . At the same time this coefficient would 
be very important for evaluating the integral in Eq. (|39|). 

There are several ways to regularize the integral in Eq. (|39| ) . Saturation effects which regularize the integral 
in the infrared could be included in the classical (multiple rescattering) approximation j4 | ^-p|, |38| , p3[ or in a 
more realistic way incorporating the effects of quantum evolution in energy 0J^,^,^| p^|j54[ "! The latter way 
involves resummation of leading logarithms of energy (i.e. powers of a s In 1/x) which arise from developing 
additional soft partons in the nuclear wave functions that will be produced in the collision. However, as 
was argued in |l^,|5^] the effect of this quantum evolution is only to produce sources of color charge off 
which the classical field would be emitted. This statement was quantified in [ p^| , p^|j5l| ] for a single nucleus. 
Unfortunately similar analysis has not been carried out for the gluon production in nuclear collisions. Thus 
we will restrict ourselves to the case of classical field emitted off the valence quarks in the nucleons of 
the colliding nuclei |35| , ^6|] . The quantum evolution will be suppressed implying that a s In 1/x < 1 and 
no extra gluons are produced. Of course generalization of Eq. (|3^) to include the full effect of saturation 
strictly speaking does not reduce to the problem of classical gluon production of J35||56| due to the non- 
Abclian nature of G a which leads to appearance of three- and four-A^ correlators in the full correlator 

^G° ly (xo)G° ly (xo)y Thus the calculations presented below should be understood as an estimate of what the 
exact answer should be in the strong field case. 
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We are going to conjecture the following procedure of generalizing the result of Eq. (J3^) to the strong field 

case of a^A^ 3 ~ 1 and o? s A\' 3 ~ 1. Let us first note that at the same order in the coupling constant as was 
employed in Eq. ( |39| ) the multiplicity of gluons produced in a collision of two quarkonium nuclei is given by 

n 

dN zMM 32 f d 2 k k x .... 

a s a o °F — / 77277 ln -T ( 44 ) 



JPbdy s 7T J (fc 2 )2 A 

with A some infrared cutoff, b the impact parameter and y rapidity of the produced gluon. As one can see 
the expression in Eq. ( f44| ) has an infrared singularity similar to the one in Eq. (|3^). Assuming that both 
singularities get regulated in a similar fashion we write with logarithmic accuracy 

<^W,(,o)>«-^^. (45) 

We thus suppose that Eq. (jig ) holds at all higher orders after inclusion of multiple rescatterings and is 
independent of whether we take a quarkonium model of the nuclei or realistic heavy ions. It has been 
conjectured in Q that the multiplicity of the produced gluons including all higher order multiple rescattering 
effects is proportional to the multiplicity of gluons in the nuclear wave function before the collision so that 
for the case of two identical colliding cylindrical nuclei 

dN c^f. (46) 



d 2 bdy a s 2n 2 

Eq. ( |46| ) has been written in p3[ ] with logarithmic accuracy. The proportionality coefficient in Eq. ( f46| ) has 
been estimated numerically to be c = 1.29±0.09 An analytical calculation of fl3£| gives c « 21n2 w 1.39 



and RHIC data seems to suggest that c = 1.23 ± 0.20 J43]. Inserting Eq. (Ug) into Eq. (45) yields 



{GIMGIM) « -c C 3 F f* 2 , (47) 



4 7T 3 a s T§ 

which, after being substituted in Eq. (|l4| ) gives the distribution of instantons 

nii(p) = no(p) cxp (- - 2 C f^ ) . (48) 

Eq. (^) is our main result. It shows that large size instantons are suppressed by the strong classical fields 
generated in the nuclear collision. As was argued in [[ll| classical fields are a good description of the produced 
gluonic medium only at the early times t Q s < 1. Thus for r ~ l/Q s the suppression starts parametrically 
when p 2 Q 2 ~ a s 1, that is our main assumption of small instantons stated in Eq. ( |l5| ) is still valid. 

Three remarks are in order here. First let us note that the instanton suppression of Eq. (^8|) results from 
Eq. (jl4|) applied to the case of nuclear collisions. Eq. flT4| ) was originally derived for the instanton distribution 
in the background of the color field generated by QCD vacuum fluctuations as considered in (2^j2^]. By 
definition of the physical vacuum the energy density in it is negative with respect to perturbative vacuum. 
Thus for the field of vacuum fluctuations (Ofy < 0. Recalling that for QCD in the chiral limit of massless 
quarks 

(K) - < G "> ( 49 ) 

we conclude that for vacuum fluctuations (G^) > 0. Therefore Eq. ( fTi] ) with the background field of vacuum 



fluctuations having positive (G^) led the authors of [^0|,^4| to conclude that the large size instantons are 
enhanced in QCD vacuum. The phenomenon is also referred to as instanton melting. The situation is very 
different for nuclear collisions. The energy density of the gluon matter produced in a heavy ion collision is 
positive with the trace of the energy-momentum tensor (Ofy > 0. With the help of Eq. (^9|) we conclude that 
produced gluon matter should have (G^) < 0, which agrees with Eq. (|47|). Thus Eq. (|l4|) with negative 
(Gf v ) leads to suppression of large size instantons in nuclear collisions as appears in Eq. (^) . 
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Secondly, it may appear that the instanton suppression of Eq. ( ]48| ) contradicts the model of soft pomerons 
in hadronic collisions proposed by the authors in fl3l[ . There each rung of the soft pomeron ladder was 
modeled by the instanton-induced transitions of two t-channel gluons into any number of s-channel gluons. If 
the instantons were suppressed so would be the effects described in ||l| . However the soft pomeron describes 
the total cross sections in the proton-proton scattering in the kinematical range of the Tevatron. There, even 
though the energies are high the colliding particles are protons (A — 1), not nuclei and unitarization effects 
are still weak if present at all. The corresponding saturation scale is presumably quite small, possibly being 
of the order of Aqcd- Our suppression of Eq. ( [48| ) would be indistinguishable from the "usual" suppression 
of large size instantons in the instanton gas models of QCD vacuum |p0| and therefore would not significantly 
alter the results of ]3l|] . Alternatively one could say that the gluon fields generated in proton-proton collisions 
at Tevatron energies are not as strong as the fields produced in nuclear collisions at similar energies and thus 
can not introduce a strong suppression of instantons. 

When the energy of the proton-proton collisions becomes extremely large saturation and unitarization 
effects begin to take place. As was argued in [0j|^,0,^3| at sufficiently high energies the saturation scale 
even in proton- proton collisions would become large, much larger than Aqcd, making most of the partons 
in the protons' wave functions perturbative. At these high energies the hard (BFKL) pomeron would be 
unitarized through multiple pomeron exchanges |B|,|9|-|l2j , which in fact would generate the large saturation 
scale Q s - At the same time unitarization of the soft pomeron is not quite understood. There is a possibility 
that the soft pomeron unitarizes through multiple pomeron exchanges, similar to the hard pomeron. Here 
we would like to outline another possibility, inspired by the soft pomeron model of [fjl) and by the result of 
Eq. ( 48 ) . The instanton-induced interactions leading to the soft pomeron behavior of the cross section would 
be suppressed in pp at high energies due to instanton suppression of Eq. (48). Thus the soft pomeron of 
pl[ could be eliminated by the large perturbative saturation scale at extremely high energies and the cross 
sections would be dominated by hard perturbative interactions which are unitary. 




FIG. 8. Distributions of instanton sizes in vacuum for QCD with three light flavors (upper curve) versus the 
distribution of instanton sizes in the saturation environment produced by a collision of two identical nuclei for c = 1 
(middle curve) and c = 2 In 2 (lower curve) with Q 2 S = 2 GeV 2 as estimated for RHIC at ^/s = 130 AGeV in |^5| . 

Third we have to note that the suppression of Eq. ( |4"§| ) is different from the instanton suppressions at 
finite T and p. There the power of the exponent is quadratic in p (see Eqs. ([}]) and (||)) while in our 
case it is quadric in p. The difference is probably due to different physical mechanisms leading to instanton 
suppression. At high temperature or density the interactions of the gluons and quarks with each other 
generate a screening mass for gluons mjj . The gluoelectric fields are therefore screened at distances roughly 
of the order of 1/mc, which leads to instanton suppression as discussed in the Introduction. In our case the 
gluon screening on the transverse scales of order 1/Q S is due to multiple rescatterings in the background 
classical field generated not by the other gluons but by the nuclei themselves. Therefore our expression in 
Eq. (|48|) is different from the high temperature and high density expressions in Eqs. (Q) and ( ft). 

In Fig. H we have plotted the perturbative distribution of instantons in vacuum no(p) of Eq. (pi) including 
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three light flavors in the beta function 6 = 9 (upper curve) together with the suppressed distributions of 
instantons in heavy ion collisions given by Eq. ( jig ) for two different values of the "liberation coefficient" c: 
the middle curve corresponds to c — 1 and the lower curve corresponds to c = 2 In 2. Saturation scale for the 
identical cylindrical nuclei was taken to be Q 2 S = 2 GeV 2 based on recent RHIC data HQ . The strong 
coupling constant at this scale was approximated by a s (Q s ) f=a 0.3. The suppression curves are plotted for 
relatively late proper time To = 1/Q S - Instanton suppression is even stronger for earlier times. 

The integral of no(p) over all p gives the total instanton density at a given space-time point. Of course 
we can not integrate the perturbative size distribution of Eq. (^) over all values of p since it diverges in the 
infrared. Instead we will employ the fit to lattice data for the instanton distribution in pure gluodynamics 
from j4(|. The density of instantons given by this fit is 10.8 fm~ . Integration of Eq. ( f48| ) with n-o(p) given 
by the fit of and with Q 2 = 2 GeV 2 , c = 2 In 2 yields us with the instanton density at t = 1/Q S of 
approximately 0.007 fm -4 . The ratio of the two numbers gives the overall suppression of instantons density 
in the central rapidity region at RHIC to be 0.0006. This estimate is made for the case of pure gluodynamics 
and the number will of course change once the effects of quarks are included. Nevertheless we still expect to 
have approximately three orders of magnitude suppression of instantons. Since the saturation scale increases 
with center of mass energy |0,^,ffl|,O| we expect the suppression to get stronger for y/s = 200 AGeV at 
RHIC and for LHC energies. Numerical estimates of [^5|,^6[ give the increase of Q s {x) for gold nuclei 
approximately by a factor of two for the LHC energy [ |15| , [16| . Therefore the expected suppression will be 
considerably stronger. 

We hope that our results open the possibility of a systematic theoretical investigation of topological effects 
in high energy nuclear collisions and introduce a different angle of looking at the problem of the interface 
between "hard" (perturbative) and "soft" (non-perturbative) interactions. 
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